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We study the quenched disordered magnetic system, which is obtained from the 2D SO (3) quan- 
tum Heisenberg model, on a square lattice, with nearest neighbors interaction, by taking a Gaussian 
random distribution of couplings centered in an antiferromagnetic coupling, J > and with a width 
AJ. Using coherent spin states we can integrate over the random variables and map the system 
onto a field theory, which is a generalization of the SO(3) nonlinear sigma model with different fla- 
vors corresponding to the replicas, coupling parameter proportional to J and having a quartic spin 
interaction proportional to the disorder ( AJ). After deriving the CP 1 version of the system, we 
perform a calculation of the free energy density in the limit of zero replicas, which fully includes the 
quantum fluctuations of the CP 1 fields 2j. We, thereby obtain the phase diagram of the system in 
terms of (T, J, AJ). This presents an ordered antiferromagnetic (AF) phase, a paramagnetic (PM) 
phase and a spin-glass (SG) phase. A critical curve separating the PM and SG phases ends at a 
quantum critical point located between the AF and SG phases, at T = 0. The Edwards- Anderson 
order parameter, as well as the magnetic susceptibilities are explicitly obtained in each of the three 
phases as a function of the three control parameters. The magnetic susceptibilities show a Curie- 
type behavior at high temperatures and exhibit a clear cusp, characteristic of the SG transition, at 
the transition line. The thermodynamic stability of the phases is investigated by a careful analysis of 
the Hessian matrix of the free energy. We show that all principal minors of the Hessian are positive 
in the limit of zero replicas, implying in particular that the SG phase is stable. 

PACS numbers: 75.50.Lk 



I. INTRODUCTION 

A spin glass is a peculiar state, which is presented by 
certain disordered magnetic materials, when the com- 
petition between opposite interactions produces what is 
known as frustration, namely the incapability of the sys- 
tem to attain the lowest possible energy state that would 
correspond to each single type of interaction. This situ- 
ation is described by ascribing to the system a random 
distribution of coupling constants, allowing interactions 
of opposite signs. As a consequence of the competition 
between different types of order, some of the spin glass 
properties are shared with paramagnetic states and some 
other with ordered, ferromagnetic or Neel states. Like 
the first ones, SG states possess zero magnetization and 
similarly to ordered magnetic states, they present break- 
down of ergodicity. The time scale of disorder is typically 
much larger than the one associated to the dynamics 
and therefore we must perform the quantum and ther- 
mal averages before averaging over the disorder, the so 
called "quenched" thermodynamical description, which 
requires the use of the replica method^. 

A SG state exhibits clear theoretical and experimental 
signatures. The former has been introduced by Edwards 
and Anderson (EA), who proposed a model for SG and 
devised an order parameter that captures one of the basic 
features of the glassy behavior, namely the occurrence of 
infinite time correlations for each spirA The same type 
of correlations exist in an ordered magnetic state such as 
the Neel state or the ferromagnetic state, for instance. In 



the SG state, however, this happens without the associ- 
ated existence of infinite spatial correlations among the 
spins and the consequent spontaneous nonzero magneti- 
zation (staggered in the case of antiferromagnetic (AF) 
order). Hence the SG state can be characterized as one 
presenting infinite time correlations (nonzero EA order 
parameter) but with zero magnetization order parame- 
ter. From the experimental point of view, one of the 
most distinctive signatures of the SG transition is a very 
sharp cusp in the magnetic susceptibility as a function of 
the temperature, right at the transition^. 

An important landmark in the study of SG was the 
derivation of a mean-field solution of a simplified ver- 
sion of the EA model, obtained by Sherrington and Kirk- 
patrick (SK)». They considered a system with classical 
Ising spins with long-range interactions, in which each 
spin would interact with any other spin in the material, 
no matter how far apart they might be. This assump- 
tion greatly facilitates, technically, the obtainment of the 
mean-field solution. Soon after, however it was realized 
that the solution of SK was unstable^. This fact has 
been generally ascribed to the so-called replica symmetry, 
which was presented by the SK solution. Indeed, later 
on, Parisi has found a stable replica symmetry breaking 
solution^. 

The long-range interactions of the SK model, however, 
are likely to be unphysical, to a large extent. Real ma- 
terials should mostly be short-range interacting, quan- 
tum SO (3) Heisenberg spin systems. Despite the large 
amount of knowledge that we have today about long- 
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range interacting SG, surprisingly, very little is known 
about the properties of short-range interacting quantum 
spin-glasses, especially with SO (3) symmetry. Apart 
from some numerical calculations^, very few analytical 
approaches exist. An appealing short-range interacting 
disordered system, which has been thoroughly investi- 
gated mostly by numerical methods is the transverse field 
Ising spin glass 8 -. Interesting related results can be found 
jjjgjio and 1 !, where a Landau-Ginzburg, phenomenologi- 
cal approach has been developed for the short-range SG 
problem. 

We have proposed a model for a disordered SO (3) 
Heisenberg-like quantum spin system with nearest neigh- 
bor interactions for which an expression for the free en- 
ergy density can be derived^. In this work, we map the 
system onto a generalized CP 1 quantum field theory in 
the continuum limit. This is a very convenient frame- 
work, because nearest neighbor interactions become just 
derivatives. At this point, however, we must be careful. 
Indeed, when taking the continuum limit of a quantum 
system Berry phases will be generated, which in gen- 
eral would not cancel when summed over the lattic o 11 ' 13 . 
We may tackle this problem by introducing disorder as 
a perturbation of an antiferromagnetic (AF) 2D Hciscn- 
berg model, for which the sum of the quantum phases 
is known to cance l 14 ' 15 . Hence, we consider a Gaussian 
random distribution of couplings centered in an AF cou- 
pling J > and with variance AJ, such that AJ <C J. 
The situation is completely different from the original EA 
model, where J = 0, and consequently the disorder can- 
not be taken as a perturbation of a Heisenberg system 3 -. 

Using the CP 1 description, we obtain a solution, which 
presents replica symmetry. Out of this, we extract the 
T x J phase diagram of the system. This exhibits a Neel 
phase at T = and p s > po, where p s = S 2 J is the 
spin stiffness (S is the spin quantum number) and po is 
a quantum critical point. This is displaced by disorder 
to the right of its original value po(0) = ^ (A = -; a: 
lattice spacing) in the pure syste m 19 ' 20 . It also contains 
a spin glass (SG) phase for temperatures below a cer- 
tain critical line T < Tc and p s < po, characterized by 
a nonzero EA order parameter and zero staggered mag- 
netization. For T > ; p s > po and also for T > Tc 
; p s < po, we find a paramagnetic (PM) phase. The 
behavior of the magnetic susceptibility, of the spin-gap 
and of the EA order parameter are analyzed in detail in 
each phase as well as on the transition line. The former 
presents a nice cusp at the transition, in agreement with 
the typical experimental behavior of spin glasses 1 . We 
also show how the phase diagram is modified by varying 
the amount of disorder 

We make a thorough investigation about the thermo- 
dynamic stability of our solution. This is done through 
a careful analysis of the Hessian matrix of the average 
free-energy density. We show in detail that all princi- 
pal minors of the Hessian matrix are strictly positive in 
the physical limit when the number of replicas reduces 
to zero. This is a necessary and sufficient condition for 



the mean-field solution to be a local minimum, hence it 
guarantees the solution is stable. Furthermore, we can 
follow the phase transition directly in the Hessian, as the 
corresponding principal minor determinants cease to be 
positive if we use the wrong solution for a given phase. 

There is an appealing physical motivation for our SG 
model, in connection to the high-Tc cuprates. Indeed, 
these materials, when undoped, are 2D Heisenberg an- 
tiferromagnets, which upon doping, develop a SG phase 
before becoming superconductors. Our model, also be- 
ing a 2D Heisenberg antiferromagnet for A = 0, describes 
precisely the AF-SG transition as we increase the disor- 
der and is therefore potentially very useful for studying 
the magnetic fluctuations of such materials. 



II. THE MODEL AND THE CONTINUUM 
LIMIT 

A. The Model 

The model consists of an SO (3) quantum Heisenberg- 
like hamiltonian, containing only nearest neighbor inter- 
actions of the spin operators Si, on the sites of a 2D 
square lattice of spacing a, 



(2.1) 



The couplings Jij are random and associated with a 
Gaussian probability distribution with variance A J and 
centered in J > 0, such that A J <C J, namely, 



1 



: exp 



(Jij - J) 2 

2(AJ) 2 



(2.2) 



We consider the quenched situation, in which, according 
to the replica method^ the average free-energy is given 
by 

F = -k B T lim - {[Z n ] av - 1) , (2.3) 
n — >o n 

where Z n is the replicated partition function for a given 
configuration of couplings Jy , 



Z n {Ji 3 } = Tr exp 



{Sf} 



a=l (ij) 



(2.4) 



and 



[Z n ] av = f ['\{dJ ij P[J ij ] Z n {J l3 }. (2.5) 
is the average thereof with the Gaussian distribution. 
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We now make use of the coherent spin states |fi(r)), 
such that 



|fi(r))=(g)(g)|nf(r)). 



i a—1 



with 



(n?(r)|S?|n?(r))=5nf(r) ) 



(2.6) 



where i: lattice sites, a: replicas, r: euclidian time, S: 
spin quantum number^, and flf (r) is a classic vector of 
unit magnitude. With the help of these coherent states 
we may express Z n as a functional integral over Slf(r), 
namely, 

Z n {Jij} = y*X>nexp j-jf%r53£ a (T)j , (2.7) 
where L a (r) 

L a (r) = L%-S 2 Y, WW ■ (r). (2.8) 



and 



(2.9) 



is the Berry phase term. 

The average over the disordered couplings Jij can then 
be performed, yielding 



[Z n ] av = J Vflexp{-Sj A } 
where Sj A is given by 



(2.10) 



S 



J, A 



(ij) 



S 4 (AJ) 2 



£ / dr^nf (r)n^(^)fi? (r)nj 6 (^^ 



(2.11) 

where summations in (a(3), as well as over the SO (3) 
components of fi, (a6) are understood. 

Using the connectivity matrix, defined as: (Ky = 1 
if (ij) are nearest neighbors and .Ky = otherwise), we 
may write the nearest neighbor sums in the last term 
of (|2.1ip as overall sums. This allows us to perform a 
Hubbard-Stratonovitch transformation that replaces the 
quartic term of (|2.11[) by 



5 A - S 4 (AJ) 2 



a, 0=1 



drdr' 



(ij) L 



(2.12) 



This is no longer a disordered system. The disorder, 
which was originally present manifests now through the 
interaction term, proportional to (AJ) 2 . In the absence 
of disorder, we would have A J — > and [Z n ] av would 
reduce to the usual coherent spin representation of the 
AF Heisenberg model, with a coupling J > 16 i 19 i 20 . 



B. Continuum Limit 

Since we are only considering the weakly disordered 
case (AJ <C J) our model, described by the effective 
action in (|2.1ip . is a perturbation of the AF 2D quan- 
tum Heisenberg model. This means we can decompose 
the classical spin f2" a (r) into antiferromagnetic and fer- 
romagnetic fluctuations as in that models. Using this, 
then it follows that the sum of the quantum Berry phases, 
Lf a , over all the lattice sites vanishes, as in the pure 
syste m 14 ' 15 . 

We can therefore take the continuum limit in the 
usual way as in the pure AF 2D quantum Heisenberg 
mode l 15 ' 16 ' 17 obtaining an SO (3) generalized relativis- 
tic nonlinear sigma model (NLSM) describing the field 
n" = (<j a , 7? Q ), which is the continuum limit of the (stag- 
gered) spin r2 Q : 



C = 



i|Vn«| 2 + ^|a T n«| 2 +zA Q (|n"| 2 -p s ) 



D r 

~2 
2 

Pt 



Q:f(r,r,r')(?:f(r,r,r') (2.13) 
<(r,T)Q^(r,T,T>£(r,r') . 



where D = 5 4 (AJ) 2 /a 2 (a: lattice parameter) and p s — 
S 2 J. The constraint n a ■ n Q = p s (no sum in a), as usual 
is implemented by the Lagrange multiplier field X a . 

This generalized NLSM contains a (A J) 2 -proportional 
trilinear interaction of n" with the Hubbard- 
Stratonovitch field (t, t'), which corresponds 
to the AJ-proportional part of (|2.11[) . namely 5a- 

Notice that a null value for J, as we have in the EA 
model^ would make the perturbation around a NLSM 
meaningless. A negative value, on the other hand, would 
correspond to the ferromagnetic Heisenberg model, which 
after taking the continuum limit, is associated to the non- 
relativistic NLSM. Here perturbation would be possible, 
however, the Berry's phases would no longer cancel. We 
emphasize, therefore, the enormous difference that ex- 
ists, both from the physical and mathematical points of 
view, in considering J as positive, negative or null in the 
Gaussian distribution of the EA model. 

The field n a = (a a ,Tr a ) is the continuum limit of the 
(staggered) spin ft a and satisfies the constraint n Q • n Q = 
p s , which is implemented by integration on A". 
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Decomposing Q a P into replica diagonal and off- 
diagonal parts, 

Q a J(r; r, t') ee 6 ab [6^ X (r, r, r') + q^(f; r, r')] (2.14) 

where q a ^ = for a = (3, we get 



J dT > [ nx 2 (T)/) + gQ/3(T)T / )ga/3(T)/) 



n Q (r) X (r,r')n Q (r') 



— n a (T)g Q/3 (T, rW) 



(2.15) 



This will be our starting point for the CP 1 formula- 
tion. In a previous worki 2 -, we took a different path. 
From (|2.15|) . we integrated over the 7?- field and thereby 
obtained an effective action for the remaining fields. 



on 

+—2" /dr'{[ X (r,r')][|zr(r)| 2 |^(r')| 5 

Ps 



where summation in i,j,a,{3 is understood. 



(3.5) 



B. The Quantum Average Free Energy 

In order to evaluate [Z n ) av in (|3.4p . we use the sta- 
tionary phase approximation. For this, we expand 
S [zf, zf*,A^, A, x(t, t'), 9 q/3 (t, t')] around the helds in 
the stationary point, in such a way that the quadratic 
fluctuations about the zf fields are taken into account, 
namely, 



III. THE CP 1 FORMULATION 

A. CP 1 Lagranagian 

We now introduce the CP 1 field in the usual way, 
namely, 



n Q (r) 



[zr(r)a ljZ f(T)} 



where the zf field satisfies the constraint 



_a|2 i |„c*|2 



(3.1) 



(3.2) 



Using the above two equations, we get the correspon- 
dence 

l|Vn«| 2 + ^|d T n*| 2 «■ 2J2 \D„z?\ 2 , (3.3) 

i=l 

where — + iA^. The above correspondence in- 
volves the functional integration over the auxiliary vector 
field Ap. 

Using (|3.ip . (|3.2[) and (|3.3|) in (|2.15[) , we may express 
the average replicated partition function in terms of the 
CP 1 fields as 



+ -J dTdr'r,r(T)M°!f(T,T')r i e(T>), (3.6) 
where rff — zf — zf s and M is the matrix 

^2 £1 ^2 £1 

S 2 S S 2 S I ' \ 6 - 7 ) 

5<°(r)^; 3 (T') 5<°(r)5^(r') 

with elements taken at the stationary fields. 

These stationary fields are such that X a (r,r) — ► A s 
and m 2 — 2iX s , which turns out to be the spin gap. Also 
X (r,r,r') - Xs (t - t') and q^(r,r,r') -> g^(r - r'). 
The staggered magnetization af is given in terms of the 
CP 1 fields as 



1 x r, a |2 , i a |21 _ 1 \ " 2 



(3.8) 



Finally, the stationary value of the gauge field is A s ^ — 0. 

Inserting (|3.6p in (|3.4p we obtain, after integrating over 
the z-fields, 

[Z n ] av = e~ nS ° a ^ "' m ' '~' A l =0 ' q ™' 5 ( T - r ').Xa(-r— (3 g) 
where 



f c S'cff k", m 2 , A s u = 0, g s Q/5 , Xsl =fivl^- [a 2 - p s ] \ 

[Z n ]av = / VzVz*VA^V X VqV\e- s , (3.4) Ls s J \2 Ls psJ J 



where S [zf , zf *, A^, A, x(r, r'), g Q/3 (r, r')] is the action 
corresponding to the lagrangian density 



-U \ P dTdr'i — 



£^ ol n ~c*|2 i -\ (\ „Ot\1 



-J, A, CP 1 



2|rJ Al zf| 2 + zA Q (|zf|^-p s ) 



+ ^ | dr' [nx 2 (r,T') + ^(r,r')^(T,r')] 



X 2 (r-r') + - g f /3 (r-r')gf /3 (r-r') 



-— [ Xs (t - r')8 a0 + qf{r - r')] <J a aA-- InDetM, 



(3.10) 
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Notice that the third term in (|3.5|) is proportional to 
n 2 (Tg and, therefore, does not contribute to (|3.10[) in the 
limit n — > 0. 

By taking the limit n — > in (|3.9p we immediately 
realize that the average free-energy is given by 

F = is cff [<, m 2 , Al = 0, g^( T - /), Xs(r - r')] . 

(3.11) 

In Appendix A we consider the determinant appear- 
ing in the last term of (|3.10p This contains the quantum 
corrections coming from the zf fields. This determinant 
runs over the i components of these fields, over the repli- 
cas and over the field configurations. We are able to 
exactly calculate the first two determinants by replac- 
ing the q a @ variables in the last term of (|3.5[) by their 
average, 



q = lim 



1 



n— >0 n(n — 1) ' 



(3.12) 



The determinant over the field configurations is most 
conveniently expressed in the space of Matsubara fre- 
quencies u> r — 2nrT, r G Z. For this purpose, we perform 
the Fourier transformation of the x's and g's. Then, using 
the expression obtained for the determinant of the quan- 
tum fluctuations (|A.10|) . we get the average free-energy 
density as the functional (we henceforth neglect the "s" 
subscript) 

/>«, to 2 , <?«>,,), xK)] - ^ [* 2 -Ps] 



-3DTY^ 



d 2 k 

(2^) 2 



In (fc 2 + M r 



Aq{uj m ) 
k 2 + M r 



(3.13) 



The EA order parameter, used to detect the SG phased, 
accordingly, is given by 

<7ea = Tq Q (3-17) 
where q~o is defined in (|3.12[) . 

C. The Stationary-Phase Equations 

By taking the variations of / with respect to the 
variables a a , to 2 , Xri f?" /3 , respectively, we obtain the 
stationary-phase equations (SPE), which are listed be- 
low. 



[to 2 - A Xo ] 5 aP - Aqf 



<t" = 



(3.18) 



ff2= ^-^E ln f 1 + ^z) +2A^g r G r (3.19) 



TA ( A 2 \ 
3DTx(-u r ) = — In 1 + — + A 2 q r G r - Aa 2 S r0 
47T V M r J 

(3.20) 



3DTq al3 {-uj r ) = A 2 q{uj r )G r + A^^-Sro (3.21) 



<J a Of3 



2n 



where 



G r — 



T 

47T 



1 



M r 



A 2 + M r 



(3.22) 



and A = 1/a is the high-momentum cutoff. 

Observe that in the absence of disorder (D, A — > 0) 
(pHD|) and (pOTj) disappear and 1(335]) . ([XT!?)) reduce 
to the well-known corresponding equations for the con- 
tinuum limit of the pure Heisenberg mode l 16 i 19 ' 20 . We 
henceforth will only consider the case fl^O. 

Equation (|3.2ip tells us that in our approximation all 
the (? Q/3 's are equal, whenever the a a, s vanish. 



where A = — and 



M r 



M(oj r ) = uj 2 + m 2 - A[xr - q r ], 



(3.14) 



where x(w r ) and q a/3 (Lu r ) are, respectively, the Fourier 
components of %(r— r') and q a/3 (r— r'). We use for these, 
the simplified notation Xr = xi^r) and q^ = q a ^(ui r )- 
From (|2.12j) . we can show that 



Qf(T,T') = (S?(T)S>:(T')) 



(3.15) 



and therefore, according to the previous decomposition of 
Q into x's and q's, we can identify xo as the static mag- 
netic susceptibility, whereas the integrated susceptibility 
is given by 



Xi = X! Xr 



(3.16) 



IV. THE PHASE DIAGRAM 
A. Paramagnetic and Spin-Glass Phases 

1. Preliminaries 

We start by searching for PM and SG phases. In both 
of them we have a = 0. Considering this fact and sum- 
ming (|3.21[) in a, /?, yields 



where T 



q(—u r )T = q(uj r )G r , 
, T ? 2 and 

7=^=3.0- 
' D \AJ 



(4.1) 



(4.2) 
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Inserting (|4~T]) in (|3~20)) . we get 

A 2 

2 2 

TO + (Jj = — r - 



or equivalently 



M r = 



A 2 



e 6np 3 (x- r -q- r ) — I 



+ A( X r-q r ), (4.3) 



(4.4) 



We can use (|4.3p in order to determine Xr — Qr- Let us 
start with r = 0. In this case, both xo and go are real 
and equation (|4.3j) is depicted in Fig. [T] 
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FIG. 1: Function appearing in the rhs of (|4.3|l (r = 0) for 
different values of the Gaussian width A J (a: 220K, b: 210K, 
c: 200K, d: 190K). The horizontal dashed line represents m 2 . 
A PM phase will only occur when this line intercepts the 
function (m 2 > mo). The physical solution corresponds to 
the left branch, m is in K and xo and qo are in K -1 



which (|4.3[) has solutions (to 2 > to 2 ,), however these will 
be such that Go < T. This is so because Go is a monoton- 
ically decreasing function of Mo and according to (|4.4|1 , 
Mo is a monotonically increasing function of to 2 , such 
that precisely Go(M (toq)) = T, as can be immediately 
inferred from (|3.22[) and from 



M (to 2 ) 



A 2 
2 



1/2 



- 1 



It immediately follows that, for to 2 > to 2 



0' 



(4.7) 



we will have 



Go(M (to 2 )) < T. Hence, the only possible solution of 



4+1). for 



> TOq , is go = 0. 



'li- 



no longer provides a solution for xoi hence we may 
now have go 7^ 0. The static susceptibility xo is now 
determined by (|4.ip . namely, Go = T. 

Now consider the case r ^ 0. We show in Appendix 
A that we always have q r — ( for r 7^ 0), because 
otherwise (|4.3p and (|4.ip again become incompatible. We 
will always have, therefore, the Xr^o determined by (|4.3p . 

Since <7ea = Tq~o it follows that go is also a SG order 
parameter and we conclude that the former phase (to 2 > 
to-o) is a paramagnetic phase (a = 0, gEA = 0), whereas 
the latter (to 2 < to 2 ,) is a SG phase (a — 0, gEA 7^ 0). 
The phase transition occurs at to 2 = TOq. 

The ratio —, which appears in the expression of the 
critical mass Too, is a measure of the amount of frustra- 
tion in the system, as we can infer from (|4.2p and the 
actual perturbation parameter. Since we are working in 
the regime of weak disorder, we take 7> 1. In the un- 
perturbed limit where the disorder is removed (7 — * 00), 
we would have TOq = and the SG phase would no longer 
exist. We see again that a disorder perturbation would 
be impossible in the original EA model, where 7 = 0. 



The function on the r.h.s. of (14.31) has a minimum at 



Xo = Xc 



6Trp s 



In 



1+7 + 



1/2 



at which the function has the value 



(4.5) 



ml 



■In 



A 2 

7 



7 



1/2 



- 1 



1 + 1 

7 



1/2 



- 1 



(4.6) 



It follows that for r = (|4.3p will only have solutions for 
to 2 > TOq. In this case, however, these solutions of (|4.3p 
are clearly not compatible with the existence of nontrivial 
solutions (g ^ 0) of P~Tj) . 

We can see this as follows. A solution go 7^ of (|4.ip 
would imply Gq = T. In the range of values of to 2 for 



2. The Paramagnetic Phase 

Let us now use the SPE in order to derive, expressions 
for the susceptibilities %i and xo in the PM phase, where 
> ml- Inserting (j3"T2%|) in ([HTT^J) . for a = and g = 



(and also g r ^o 
susceptibility 



0), we readily obtain, for the integrated 



(4.8) 



It follows from (|3~20)) that 



X™ = ^ T— E t ln ( A2 + M r) - kM P ]. (4.9) 

The previous sums are dominated by large values of 
uj r . In this case, we show in Appendix A that (|4.3p or 
equivalently IjB.ip and (|B.2p yield the solution 



X(Vr) 



A 2 



6irp s (to 2 + w 2 ) 



(4.10) 
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We may then evaluate the two sums above, obtaining for 
them an explicit expression 

HA 2 + M r) ~ lnM r ] = lnT(m 2 ,T, 7 ), (4.11) 



where 



sinh 



2T 



sinh" 



/-X- 
2T 



sinh 2 (^^^Jm 2 + sinh 2 



,2 _ Af 



v/7 



™ 2 + 7f K"^ 



(4.12) 



with 

X± = (A 2 + 2m 2 



1 1 

T 



A 4 



1 



2 1 2y (A 2 + 2m 2 ) 2 

(4.13) 

The static susceptibility in the PM phase, therefore, is 
given by 



i 



i 



In T(m 2 , T, 7). 



3T 67rp s 

The function T(to 2 ,T, 7) has the properties 



T(to 2 ,T, 7 ) 



T»A 



and 



In T(m 2 , T, 7) 

where (for 7 >> 1) 
A 



T^O 27T 



r 



PO, 



We also have 



1 + 1 

7 



l + iln(l+ 7 ) 



— In T (to 2 , T, 7) 



Xc 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



where x& is given by (|4.5| . implying that the critical 
value of x™ is Xcr- 

We see that Xo satisfies the Curie law at high- 
temperatures 



pm t»a J_ 

XO Qrp 



and diverges as 



PM T-+0 J_ 
AO 3X" 



Pa. 



(4.19) 



(4.20) 



for T — > 0. As we will see this is the expected behavior 
for (p s > po)i where an AF phase appears at T = 0. For 



(p s < po), conversely, we will see that the PM-SG phase 
transition occurs at a finite T c (Fig. [5] and Fig. [3]) and 
the previous expression is no longer valid. 

Notice that the following general relation involving 
the integrated susceptibility^ is automatically satisfied 
by xf M and q PM (uJ r ) = 0: " 



Xi = 



(4.21) 



5. Tfte Spin-Glass Phase 

We now turn to the the SG phase (to 2 < m§). Now 
4.ip . or equivalently Go = T implies 



M 



A 2 
2 



1 + 



7 



1/2 





_ A 2 






1 




1- 






7 







which coincides with (|4.7[) . 

Then (Prgj) yields (for 7 » 1) 



'in 



1 1 - „, 2 _ . 1 



SG - i^~7T— lnT(TO^,T, 7 ) 



(4.22) 



Tn(l + 7). (4-23) 



3T 67rp s 67rp. 
From this and (|4.22p we obtain 

Xn G = — In T (m 2 , T, 7) — — [to 2 - to 2 ] . (4.24) 

Ao 3T 67rp s V ' ' " 2D 1 J V ' 

Since (pHO)) still holds for r ^ 0, we still have 



52 X(^r) = 7— lnT(TO 2 ,T, 7 ) 
J^o 6?r ^ 



and therefore we get 



1 



V SG - — - -^fTO 2 - TO 2 1 

Xl ~3T 2D [m ° mJ 



(4.25) 



(4.26) 



Comparing ([424)1 and (|4~T4)) and also ((426)) and ((48)) 
we can identify a clear cusp at the transition appearing in 
both susceptibilities. This is an important result, since 
the presence of these cusps is a benchmark of the SG 
transition and has been experimentally observed in many 
materials presenting a SG phase 1 . 



B. The Neel Phase 

Let us now search for an ordered Neel phase, for which 
(7^0. We see that in this case the quantity between 
brackets in (|3.18p must vanish. By summing it in a, (3, 
we conclude that Mo = in this phase. According to 
p,19p . however, this can only happen at T = 0, oth- 
erwise we would have an unphysical infinite imaginary 
staggered magnetization a. This is in agreement with 
the Mermin- Wagner theorem—, and is a clear evidence 
that our approach goes beyond mean-field. 
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Now, Mq = implies 



Critical Curve 



% 



m 



(4.27) 



On the other hand, for a ^ 0, (pH!?)) and (|3~2T))) imply, 
instead of 



xf F 



2a^_ 

Ps 



From this, it follows that 



v AF - — 
Xo - 3T 



1 



6-Kps 

which in the limit T — 



af t^o J_ 

Xo ~ 3T 



lnT(m 2 ,T, 7 ) 
reduces to 
1 



2q 2 
3Tp s 



Po 

Ps 



2a 2 
3Tp s 



(4.28) 



(4.29) 



(4.30) 



for p s > p Q . 

On the other hand, (|4.2ip leads to 



o AF - - 
% —j, 



2a 2 

Ps 



(4.31) 



Equations (|4T27|) (j4730|) and (|4~3"Tjl . allow us to solve for 



,.AF ^AF 
Xo 



and a. We get, for T 
, 1 



o. 

[Ps ~ Po] 



and 



4 F 



i 



[Ps - Po] 



(4.32) 



(4.33) 



This implies, according to (|4.27p a zero spin gap: m 2 



0. 



From l(4~28| and (|4~32"1) we obtain 



Xi 



AF 



1 

12T 



Po 

Ps 



(4.34) 



Notice that the susceptibilities diverge for T — ► as 
they should. The EA parameter, however, remains finite: 



n AF 

Iea 



Po\ 



(4.35) 



We see that, indeed, there is an AF phase characterized 
by ( a 7^ 0j Iea ^ 0) on the line (T = 0,p s > po), with 
po given by (|4.17p . In the absence of disorder (A J = 
0,7 — > oo), po —* Po(0) = which is the well-known 
quantum critical coupling determining the boundary of 
the AF phase in the pure 2D AF Heisenberg model at 
T = Oi&i^. The effect of disorder on the AF phase is 
to displace the quantum critical point (QCP) to the right. 
This result should be expected on physical grounds: in 
the presence of disorder a larger coupling is required, to 
stabilize an ordered AF phase. 



We have seen that the parameter m 2 determines the 
transition between the PM and SG phases. It is impor- 
tant, consequently, to see how it depends on the con- 
trol parameters of our system, namely, T, J and A J, or, 
equivalently, T, p s and 7. We may obtain an equation 
for to 2 by using (|4~2l"1) . (|4~23]) . (142%]) and the fact that 
Qrjto = 0- These yield 



2D 



m — ms = 



Xc 



3p s [ 3T 6irp 
From this we get 



^ + 7T^lnT(m 2 ,r, 7 ) 



SG 



= *£l r. 

2D L 



(4.36) 



(4.37) 



for m? < mg. We see that % —* at the transition as 
it should. For m 2 > to 2 , we have M = as seen above. 

We may determine the critical curve by observing that 
the critical condition to 2 = mg implies 



1 

37; 



671-/0. 



lnT(mo,T c ,7) = Xc 



(4.38) 



For T c <C A, which corresponds the situation found in 
realistic systems, this becomes, near the quantum critical 
point (p a < po), 



T 

- 1 c 

2^ 



In (A) lu(l + :7 



Po 



(4.39) 



which is the equation for the critical curve separating the 
PM and SG phases. Notice that it meets the T = axis, 
precisely at the quantum critical point po, separating the 
SG from the AF phase. We plot the T c x p s phase diagram 
corresponding to (|4.39|) . for a fixed value of 7 in Fig|2] 

In Fig. [31 we plot again T c x p s , for a fixed value of 
the amount of disorder, namely, the Gaussian width A J. 

Now we plot in Fig. 2] the associated T c x A J phase 
diagram, for different values of the spin stiffness. The 
latter shows the AF-SG transition as a function of in- 
creasing disorder. This is the type of transition which is 
observed in the high-Tc cuprates. In order to describe 
it, we must relate the doping parameter of these mate- 
rials to our disorder parameter, AJ. We are presently 
investigating this point. 

Below, we plot the quantum critical point po as a func- 
tion of the amount of disorder, for different values of p s . 



D. Critical Behavior 

The critical behavior of relevant quantities may be 
determined by analyzing the function Y(to 2 ,T, 7) for 



9 




0.96 0.98 



Ps / P0 (0) 



20 



15 





PM 


(a)/ 

/ (b y 
/ / (c) / 


AF 




/ / SG 



100 



200 



300 



400 500 

A.T 



600 



700 



800 



FIG. 2: Phase diagram for a fixed value of 7. (7 = 10 2 , A = 
10 3 ). The critical curve corresponds to (|4.39|l and is valid near 
the QCP p (solid curve). p (0) = A/2vr is the QCP of the 
pure AF system. Notice that disorder besides creating the SG 
phase, displaces the QCP to the right. The value ascribed to 
A is a realistic one in K (A — > '-^ A; v s : spin- wave velocity). 
The resulting temperatures naturally appear with the correct 
order of magnitude, in K, found in real SG systems^. 
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FIG. 3: Phase diagram for A = 10 3 and AJ (a: 200K, b: 
350K, c: 460K). Notice that the AF phase (thick line) is dis- 
placed to the right as we increase the disorder. 



T r~j T c and m 2 ~ m§. This yields, near the transition, 
for p s < po, 



^-e^: lnT(m2 ' T ' 7) ^(|) v 



and 



AttA 



T-T c 



[P0 - Ps 



(4.40) 



(4.41) 



FIG. 4: Phase diagram T c x AJ, for A = 10 3 and p s /po(0) (a: 
1.01, b: 1.03, c: 1.05). Notice that as we increase the value of 
the center of the Gaussian (p B ), a larger amount of disorder 
will be required for the SG phase. 




200 



FIG. 5: The quantum critical point as a function of the 
amount of disorder for different values of the average cou- 
pling, p s ( a: 180K , b: 200K, c: 220K ). Notice that in the 
limit of zero disorder po — * Po(0) = A/2-7T. 



the SG order parameter and susceptibilities for T > T c 
and p s < po : 



and for T <T C and p s < p : 



D 



T c -T 



[ Po - p.] (4.43) 



From these expressions and (|4~TI)) . Igjgj) , (|g^3|l . (|4^5|) 
and (|4.23p we can fully determine the critical behavior of 



sg _i 2?r Ap s 

Xl ~ 3T D 



T n — T 



[po 



(4.44) 



10 



-SG 
% 



6irAp s 
D 



T c -T 



[PO - Pa], 



(4.45) 



follows. 

We plot Xo m Fig. [U We can see the characteristic 
cusps of the SG transition occurring in these magnetic 
susceptibilities. 
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F. Effect of Quantum Fluctuations 

As mentioned before, the effect of quantum fluctua- 
tions, introduced by integration over the quadratic fluc- 
tuations of the CP 1 fields and embodied in the last term 
of (|3.13[) is essential for the obtainment of a sensible 
solution for this system. Should we do a pure mean- 
field approach, by disregarding these fluctuations, we 
would obtain a Neel state for any temperature, in ob- 
vious disagreement with the Mermin- Wagner theorem. 
Conversely, taking these fluctuations into account washes 
out the Neel phase to T — 0, but leaving a spin glass state 
below the critical curve. The fact that the AF phase is 
removed to T=0 shows that our approach transcends the 
mean-field approximation. The SG phase is robust to 
such fluctuations. In the spirit of the loop expansion that 
is being done, we assume that higher quantum fluctua- 
tions will not change this picture qualitatively. Anyway, 
considering only the quadratic fluctuations is always a 
valid first approach to a difficult problem. 



THE THERMODYNAMIC STABILITY 



FIG. 6: The static susceptibility for different values of p s : 
p s /po(0) =(a)1.005 ,(b)0.993, (c)0.981. The cusps, character- 
istic of the SG transition, occur at the corresponding critical 
temperatures (in K). \o is in K -1 . 

For p s > p and T > 0, we always have m 2 > m 2 , and 
& = q~o = i.e. the system is in the PM phase for any 
finite temperature. 



E. Dependence on Disorder 

Let us now examine the A J dependence of the phase 
diagram. For this we make AJ — > AJ(1 + e), with |e| <C 
1, for fixed p s , and study how the relevant quantities 
change. We find, 



For p s < po and a fixed p s , we also obtain 
Tf - T r . \mlY 



(4.46) 



inApo 



°-^ P °- £ -^ln 7 ., 



Ap 



2iT~/Ap 



where Ap = po — p s . 

The SG order parameter changes as 



(4.47) 



(4.48) 



We see that increasing the amount of disorder, through 
an increment of the Gaussian width, will increase po and, 
for a fixed p s , also T c , m 2 , and qg G . Conversely, decreas- 
ing the amount of disorder by narrowing the Gaussian 
width will produce the opposite effects. 



We finally consider the important question of the ther- 
modynamic stability of the phases. As is well known the 
replica symmetric mean-field solution of the SK model 
turned out to be unstable. It is, therefore, absolutely 
necessary to analyze the stability of any solution to a 
SG. We will focus on the SG and PM phases. The sta- 
bility of the AF phase should not be a problem. 

For studying the stability of the solution, we must con- 
sider the Hessian matrix of the free-energy density 



/ = f[a a ;q^(uj ), ...q afi (u r )...; A; X M, xK), ■••], 

(5.1) 

where the variables q a °(u) r ) and x( w r) are complex, such 

that X-r = Xr an d Q-r = 1 
The Hessian is given by 



a/3 



a 2 / 



(5.2) 



where fa = a a , q a P (uj ) , Re q ap (w P ) , Im q^ (w P ) , A, x(w ) , 
Re x(aj r ),Im x(w r ), where the index r runs from 1 to 
oo. This is a matrix with entries of dimensions 
lim„^ [n;n(n - 1) , ...,2ra(n - l) r , 1; 1 , -, 2 r , ...] 
corresponding, respectively, to derivatives with respect 
to each of the above variables. Since we are interested 
here in the SG and PM phases, we will take (|5.2p at 

(7 = 0. 

The elements of the Hessian matrix are as follows: 
there are four overall diagonal and six crossed elements, 
namely aa, qq, AA, XX, eg, aX, ax, <?A, qx and A%- 

The a a term is, 



d 2 f 



1 



da a daP 



~M S ai:s - Aq C afj ] 



(5.3) 
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The q r q s terms are proportional to S rs . For r = 
d 2 f 



dqfdq^ 
where 



[aoS^S^ + 6 C" 7 C^] = K , (5.4) 



ci r — A. — . b r — 



[n{n-l)f 



and 



(5.5) 



where F r = G r + H r 
For r ^ 0, we have 



d 2 f 



dRc XrdRe \s 



d 2 f 



aim Xrdlm Xs 



A 2 [T + Re F r } = L r (U)S rs , (5.12) 



^l 2 [r - Re F r ] = L r (22)5 rs , (5.13) 



T _ 

H r = A—q r 
Air 



1 



1 



M 2 (A 2 + M r ) 2 
For the r ^ blocks, we have four terms, namely 



d 2 f 



(5.6) dRc xrdlm Xs 



A 2 [lm F r ] = L r (l2)S rs = L r (21)6 r 



(5.14) 



The aq, a A and ax crossed terms vanish for a = 0: 



d 2 f 



dRc q a 0(u r )dRe q a ^{oj s ) 



[arS^S^ + Re bX^C^] = K r (ll)S rs , (5.7) 



d 2 f 



d 2 f d 2 f 



da<*dqf da a 0\ da a d Xr 



= 



(5.15) 



The only non-vanishing qX and qx terms are the ones 
for which r — 0. These are respectively, 



d 2 f 



dim q a P(uj r )dIm ««0(w 8 ) 



[a r 6 a ~«5^ - Rc brC^Cf"] = K r (22)S rs , (5.8) 



d 2 f 



dRc q a P(u r )dIm q al3 {io s ) 



[-Im b r C cn C f3r '] = K r (12)S rs = K r (21)5 rs . (5.9) 
The AA term is, 



d 2 f = 

<9A<9A dm 2 dm 2 



d 2 f 



Yp r + H r ] = y (5.10) 



The XrXs terms are also proportional to S rs . For r = 0, 
we have 



d 2 f 
dxodxo 



A 2 [r-F ] = L 



(5.11) 



d 2 / . d 2 f 

= i - jL - - = -iAC 



dXdqf dm 2 dqf 



d 2 f 



where 



dxadqo 



Co 



a/3 



-A 2 C 



H 



n(n — 1) 



(5.16) 



(5.17) 



(5.18 



The r ^ terms vanish because q r = (and consequently 
H r = 0) for r ^ 0. 

Finally, the Ax terms are 



5 2 / 



a 2 / 



<9A<9Re Xr dm 2 dRe Xr 



-iA[Re F r ] (5.19) 



and 



d 2 f 



d 2 f 



9A<9Im Xr dm 2 dim Xr 
The complete Hessian is as follows 



iA[lmF r }. (5.20) 
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Cq/3 : 







•■■ 



K 



-iAC 
-A 2 C 



K r 

-iAC Q ■■■ ■•■ 
-A 2 C 

••• 
••• 



-iAC -A 2 C 



-iAC -A 2 C 





-iAF 

-iARe G r 
-iAIra G r 











-iAF ■ ■ ■ -i^4Re G r —iAlm G r 



L r (ll) 
L r (21) 



L r (12) 
L r {22) 



r 



Here if r is the matrix block with elements given by 



A necessary and sufficient condition for the mean field 
solution to be a local minimum is to have all the principal 
minors of the Hessian positive. This would be equivalent 
to having all the eigenvalues of the Hessian positive. The 
principal minors are the determinants of the matrices ob- 
tained from the original matrix by successively striping 
the last line and the last column, starting from the ma- 
trix itself and ultimately reaching the (11) element. In 
the n — ► limit we have the following set of principal mi- 
nors of H: D(j , D qo ,..., D qr ,.,., D\, D Xo ,..., D x * r , D x »... . 
In the previous expressions, the prime and double prime 
refer to the two principal minor determinants generated 
by the sub-matrix L r (ij). 

We have carefully evaluated each of these determinants 
(see Appendix C) in the limit n — > 0, for a — (PM and 
SG phases). 

Defining 



\G S 



with 



and 



P = lim P r 

r — >oo 



(5.22) 



(5.23) 



(In Appendix D, we demonstrate that both P and G 
are finite and positive), 
we obtain 



D a = 1; 



D xa = A 2 



G()Dq + G 



(r - Go) 



2 1 



P: 



D' Xr = A 2 [T + Re G r ]D' Xr i - 



A 4 



r 

and 



(r-G ) 2 [Re G r } 2 



l[A 2 [T 2 -\G S \ 2 ] 



(5.25) 



(5.26) 



D qr =D qo P r ; D x = D qo PG; (5.27) 



(5.28) 



P (5.29) 



Dl = [A^T 2 - \G r \ 2 }} D'l_ x 
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^(r-G ) 2 A 2 |G r | 2 [r-ReG r ] 



or equivalently, 



Y[a 2 [r 2 -|G s | 2 ] 



p 



(5.30) 



^=n^ 4 [ r2 -i^i 2 ]]^xo- 



.4 2 



— (r-G ) 2 ^A 4 |G s | 2 [r-ReG s ] 



s=l 



[r 2 -|G t | 2 ] 

(5.31) 

We see that all the principal minors, except D a = 1 
can be written in the form 



P 



£D qo + r)[T - G f 



(5.32) 



where £ is a positive factor and 77 is either positive or zero. 
These properties follow from the fact that T > \G r \, for 
r 7^ 0, as we demonstrate in Appendix B. 

From (|5.32[) , we see that D qo is a key piece in the eval- 
uation of the principal minors. In the PM phase we have 
<7o = 0, Hq = and Go < T, therefore it follows that 
D qo > 0, hence all the principal minors are positive. In 
the SG phase, conversely, q~o > 0, Hq > and Go = T, 
implying that again D qo > 0. It follows that all the prin- 
cipal minors are positive also in the SG phase. At the 
transition, all principal minors vanish, except for D a = 1. 

The above result establishes the thermodynamic sta- 
bility of the SG and PM phases obtained from our solu- 
tion. Furthermore, we can see the phase transition oc- 
curring by direct inspection of the principal minors, since 
they vanish at the transition point. This result rules out 
the usual instabilities found in the long-range interact- 
ing replica-symmetric solutions. We can still have, of 
course, meta-stability, which seems to be a characteris- 
tics of spin-glasses. This would deserve further investiga- 
tion. For this purpose, a promising procedure would be 
to use the method of quantum annealing 21 , in order to 
find the absolute minimum of the free energy. 



VI. CONCLUDING REMARKS 



we would obtain a Chern-Simons term for the ^4^-field in 
the final CP 1 version of the model. This system deserves 
a deeper investigation but, presumably, the presence of 
this term will not modify the phase structure found here. 

Our solution takes into account the quantum fluctua- 
tions of the CP 1 fields and therefore transcends the mean- 
field approximation. This fact becomes evident, when we 
note that our solution does not predict any ordered AF 
phase at T ^ 0, in agreement with the Mermin- Wagner 
theorem, but contrary to what a mean field approxima- 
tion would yield. 

The stable SG phase derives from a replica-symmetric 
solution. This indicates that, in the case of short-ranged 
interactions, there is no basic clash between the replica- 
symmetry of the SG solution and its stability. We are 
naturally led to inquire, therefore, whether the instabil- 
ity, which has been found in the SK solution, is actually 
produced by the long-range interaction itself, rather than 
by the replica symmetry it possesses. 

The plots of the magnetic susceptibilities versus the 
temperature, exhibit the characteristic cusps, experimen- 
tally found at the PM-SG transition, in materials exhibit- 
ing the SG phase. By choosing a realistic value for the 
momentum cutoff A, we see that the cusps occur precisely 
at the temperature values, which are observed experi- 
mentally. This provides clear evidence that our model is 
really capable of describing realistic SG systems and our 
results are not an artifact of the mean-field approxima- 
tion. 

Our model nicely describes the AF-SG transition, 
which occurs as we increase the amount of disorder, hence 
it will be probably useful in the description of the cor- 
responding transition in the high-Tc cuprates. We are 
currently investigating this point. 
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We have proposed a model for describing a short-range 
interacting, disordered quantum magnetic systems with 
SO (3) symmetry on a square lattice. The random dis- 
tribution of couplings is a Gaussian biased to an AF 
coupling. A replica symmetric stable solution was ob- 
tained, which clearly shows the existence of a stable gen- 
uine SG thermodynamical phase, at a finite T. This can 
be seen directly from the solution, but also by examin- 
ing the principal minors of the Hessian matrix of the free 
energy, which are all positive, both in the PM and in the 
SG phases, but vanish at the transition. 

The use of AJ < J > 0, allowed us to assume the 
cancelation of the Berry phases. Relaxing this condition, 



APPENDIX A: DETERMINANT OF THE 
QUANTUM FLUCTUATIONS 

Let us evaluate here the determinant of the matrix of 
quantum fluctuations M appearing in (|3.10p . From (|3.7[) 
and (|3.5p we obtain 



K°P + A 



a/3 
11 



n af3 
°12 



c, 



21 



K a P + A 



a0 
22 



(A.I) 
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where, already in momentum-frequency space, 



K 



a/3 _ 



\k\ 2 + lu 2 + m 2 



S a0 



(A.2) 



APPENDIX B: xK) AND q(u r ) FOR r ^ 
1. The Stationary Phase Equations 



and 



ACtP _ Ad/3 * 

ij ~ 71 Z i Z J> 



C?f = A^z iZj . 



In these expressions, 



A af3 



2D 



(A.3) 



(A.4) 



[[q(u r ) - xM]S aP - q(u r ) C^] , (A.5) 



where C a/3 is the n x n matrix with all elements equal 
to one and x(w r ) and q(u r ) are, respectively, the Fourier 
transforms of %(r — r') and q(r — r'). 
We want to calculate 



Let us define 67rp s [x((xv) — <z(k> r )] = a + id, for an 
arbitrary r ^ 0. Then, the mean field equation (|4.3|) 
yields the two equations 



29 [cosh a — cos 9] — 7 sin 0. 



and 



m + uj,. 



A 2 [cos ( 



A 2 



2 [cosh a — cos 9] 



-a. 



(B.l) 



(B.2) 



These equations imply |G r | < T, as we demonstrate 
below. 



In Det M = In Det det det 1 



(A.6) 



where the three determinants run, respectively over the 
momentum- frequency arguments of the fields, the repli- 
cas and the z-field components. The first determinant 
can be easily evaluated by diagonalizing (|A. 1[) . We get 

In Det M = In Det det \K ap + p s A a0 ] (A.7) 
fc Wr <*P 1 

The two remaining determinants were evaluated irJ^.. 
The one over the replicas yields exactly, for n — > 0, 



det \K af> + Ps A afS ] = Nq 



1 Nl ' 

1 — n — 

N . 



(A. 



where 



and 



N = |k| 2 +^ 2 + m 2 



2D 



(Xr 



Ni 



2D 



Inserting (|A.8|) in (|A.7|) . finally, we can write the last 
determinant as trace, which for n — > reads 



In Det M = n Tr 

kuj r 



Ni 

In N 

N 



This is given by 



In Det M = iiFV 



d 2 k 
(2nf 



In (k 2 + M r ) 



where A= 2 ^- and M r is given by (|3.14[) . 



(A.9) 



Aq{u m ) ' 
k 2 + M r _ 
(A.10) 



2. |G,.|<r 



Indeed, from (|4.3p . we have 
A 2 



This yields, from (|3~22"j) 



2r 



e a+i6 _ l 



G r = — [cosh(a -I- i9) - 1] 

7 



which implies 



G r | 2 = - [2 (cosh a- cos 6»)f 



From (|B.ip . we see that, for 6* 7^ 



IGU =r 



< r 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



For 9 — 0, conversely, G r is real and we see from Fig.l 
that a r < ao (for r ^ 0). Since G r is a monotonically 
increasing function of a ri it follows that G r < Go- Now, 
as we have seen, Go < T, hence G r < T for = 0. 
We conclude therefore that, for r ^ 0, we always have 

|G r | < r. 

For a non- vanishing g r , however, Eq. (|4.ip implies 
|G r | = T. We conclude, therefore, that (|4.3p . or equiva- 
lently (|B. 1[) and (|B.2[) will only admit solutions for r ^ 
when q r — 0. In this work, therefore we always have 
q r ^o = 0. We also choose the 9 = solutions of (|B.1[) 
and (|B.2[) for all r ^ 0, which imply real Xr's. 
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APPENDIX C: THE PRINCIPAL MINORS OF 
THE HESSIAN 

1. The Determinant D a 

According to (|5 . 3[) . the determinant of the (era) n x n 
block of the Hessian, for finite n is given by 



D a (n) = 

n 



n n 



In the limit n — > 0, this gives 



90 



D a = lim = 1 - nA— = 1, 



(C.l) 



(C.2) 



where we used the fact that 



n— >o \ n 



The Determinant _D„ 



From (|5.27p ). we have the determinant of the (go<Zo) 
n(n — 1) x n(n — 1) block of the Hessian, for finite n, 
given by 



D qo (n)=a Q n{n 1} 11 [a - n(n - l)b ] , 

where clq and bo are given by (|5.5p , for r = 0. 
In the n — > limit, this gives 



(C.3) 



4. The Determinant Da 
From (|5.2[) . we have that 

oo 

D x (n) = ^+A 2 n(n-l)C 2 4" (n - 1) - 11 JJ 

r^O 

where we used the fact that c r = for r ^ 0. 
Now, using (|5.5jl and (|5.18|) . we get 



r 2 - |G r 



lim A 2 n(n - l)Cg4 n(n - 1)_11 



r 2 

(C.7) 



(C.8) 



Inserting in (|C.7[) and using ()5.10j) . we immediately ob- 
tain D\, considering that always Hq[T — Gq] = 0. 



5. The Determinant D xa 



From (|5.2p . we get 

£> Xo (n) = A 2 [r - (Go + Hq)] D x + A 2 F Q 2 D qoa 

+ A i {2F Q - <p)n{n - l)C 2 a n(n_1)_11 P (C.9) 



Taking the limit n — > 0, using (|C.8[) and (|5 . 1 0[) . we obtain 
(|5.28p . after a little algebra. 



6. The Determinants D' X1 and D' X1 



D qo = 



lim (1) 
which is (|5.26p . 



T Go — H 
n n(n — 1) 



= ~[(r-G ) 



Ho], 
(C.4) 



From (|5.2p . using the fact that H r = for r ^ 0, we 



get 



L> Xl (n) = A 2 [r + ReGi]^ 



The Determinant D 9r ,r / 1 



+,4 2 (Re Gi) 2 [,4 2 [r - Go - H ]D q 



From (|5.7p - (l5.9p , we can show that the determinant of 
the 2n(n — 1) x 2n(n — l)-dimensional, (q r q r ) block of the 
Hessian, for finite n, given by 



D Q r \n) = U 2 ^"" 1 )] + 2n(n-l)& 4 2n(n ~ 13 ~ 1] W 



where a r and 6 r are given by (|5.5p . 

The limit n — > can be taken in the same way as 
we did in the previous subsection. Using the fact that 
H r = for r ^ 0, we obtain, 



r 2 - \G r 

r 2 



(C.6) 



Since the g-part of the Hessian is block-diagonal, we 
immediately establish (|5.27p for D qr . The limit r — > oo 
exists, as we show in Appendix C. 



-A 4 2n(n - l)Cg4 M "" 1)_11 



P. 



(CIO) 



Using (|C.8p . after some algebra we obtain, in the limit 
0, 



; [r-G ] 5 



D' xi = A 2 [T + Re Gi] D Xo + A 4 (Re Gi) : 



Following an analogous procedure, we obtain 



(C.ll) 



d>; i= a* [r 2 -\ Gl \ 2 ] d xo + 



A 4 |Gi| 2 [r-Re Gy 



[T - Go] 2 



(C.12) 



Considering the cases of -D X2 and D x , , it is not difficult 
to obtain, by induction, the general expression for D' x 
and . eas. (lQ9>((OT] l. 
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APPENDIX D: FINITENESS OF G AND P which is finite. It follows that < G < oo. 



Theorem: < G < oo 



From (J4~4f and (pT22|) . or directly from (|B.5[) . after 
choosing the solution # = 0, we can write 



G r 



T 



2ttA 2 



[cosh(67rp s Xr) - 1] > (D.l) 



As argued before, according to (|4. 3|) or (|4.4p . for large 
values of r, we have \r given by (|4.10p . In this case, 



G r 



T 



4ttA 2 



A 2 



l 2 



m 2 + uj- 



(D.2) 



and, for a sufficiently large but finite N , G can be written 
as 



JY 



G — Gq -\- 2 ^ ^ G r 



T 



r=l 



2ttA< 



E 

r>N 



A 2 



TO 2 + W 2 



(D.3) 



The first two terms are obviously finite and positive. The 
third term, is clearly smaller than 



T 



4ttA 2 



E 

r=0 



A 2 



m 2 + w 2 



(D.4) 



2. Theorem: < P < 1 
From (j5T2"2"]) and (folio) , we have 

OO 

hiP = J^ln 



> -4 



(D.5) 



Using the same idea of the previous subsection, we can 
write, for a sufficiently large but finite N 



N 



InP = J^ln 



1-1^ 



r>N 



A 2 



to 2 + cj 2 



(D.6) 

The first term is obviously finite and negative. The mod- 
ulus of the second term is clearly smaller than 



^2 E 

' t — n 



A 2 



m 2 + cj 2 



(D.7) 



which is finite. It follows that < |lnP| < oo, with 
InP < 0. We conclude, therefore, that < P < 1. 
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